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Abstract 

Perturbative a' corrections to the low energy string effective action have been recently found 
to have a potentially regularizing effect on the singularity of the lowest order pre-big-bang 
solutions. Whether they actually regularize it, however, cannot be determined working at 
any finite order in a perturbative expansion in powers of the string constant a', because 
of scheme dependence ambiguities. Physically, these corrections are dominated by the in- 
tegration over the first few massive string states. Very massive string modes, instead, can 
have a regularizing effect which is non-perturbative in a' and which basically comes from 
the fact that in a gravitational field with Hubble constant H they are produced with an 
effective Hawking temperature T = H/{27c), and an infinite production rate would occur if 
this temperature exceeded the Hagedorn temperature. We discuss technical and conceptual 
difficulties of this non-perturbative regularization mechanism. 



CERN-TH/97-228 
August 1997 
revised April 1998 



^Permanent address: INFN and Dipartimento di Fisica, piazza Torricelli 2, 1-56100 Pisa, Italy. 



1 Introduction 



String theory is an appropriate framework for discussing the singularities of general relativity 
and in particular the big-bang singularity, and a large number of works have been devoted 
to the study of cosmological solutions of the low energy effective action of string theory (see 
e.g. [1-5]). At lowest order in the string constant a', the solutions of the equations of motion 
still reach a singularity at a finite value of time, say t = 0, when we evolve the present state 
of the Universe backward in time 0. This is not surprising since, in the lowest order effective 
action, a' only enters as an overall constant and therefore drops from the equations of motion; 
therefore, there is no scale at which the singularity can be regularized. For homogeneous 
fields, the low energy action has a symmetry, scale factor duality, which relates different 
solutions of the equations of motions [|l], In the simplest case of an isotropic Friedmann- 
Robert son- Walker (FRW) metric with scale factor a{t), and of a vanishing antisymmetric 
tensor field B^^, it reads 



where is the dilaton field and d is the number of spatial dimensions. This symmetry can be 
generalized to a global 0{d,d) symmetry in the more general case of non-diagonal metrics 
and non- vanishing B^,^. Since it involves the dilaton field, it is a string symmetry, with 
no counterpart in Einstein gravity. Combining scale factor duality with time reversal, it is 
possible to associate with every 'post-big-bang' solution defined for t > a 'pre-big-bang' 
solution defined for t <0, through the transformation a{t) l/a(— t). Of course, at lowest 
order in a', the pre-big-bang solution also runs into a singularity as t ^ 0^. It is therefore 
natural to ask whether the inclusion of corrections allows a smooth transition between the 
pre-big-bang and the post-big-bang solutions, providing a non-singular cosmological model. 

The effective action of string theory has two different types of perturbative corrections: 
higher orders in a', which are genuinely stringy corrections related to the finite extent of the 
string, and loop corrections, which carry higher powers of e"^. Perturbative a' corrections 
provide a scale for the regularization of the singularity and have been studied in 0] , where 
it has been found that the equations of motion in the case of constant curvature and linear 
dilaton reduce, at all orders in a', to a set of {d+ 1) algebraic equations in {d+ 1) unknowns. 
If these algebraic equations have a real solution, this can act as late time attractor of the 
lowest order pre-big-bang solution, as has been shown on specific examples at 0{a'). The 
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singularity is then replaced by a phase of De Sitter inflation with linearly growing dilaton. 
To complete the transition to the standard radiation-dominated FRW model (the graceful 
exit problem) O^e'^) corrections must play a crucial role [|, ||, |1^, [Tl[ . 



In this paper we examine again the role of a' corrections in string cosmology. In sect. 2, 
after recalling the results of ref. [0] and extending them to the case of compact dimensions 
(sect. 2.1), we discuss in some detail the scheme dependence of the results (sect. 2.2) and the 
relation with scale factor duality (sect. 2.3). Unfortunately, to verify that the mechanism 
proposed in [0] does take place would require the knowledge of a beta function at all orders 
in a'. Therefore, in sect. 3 we try to obtain a better understanding of the physics behind 
the perturbative corrections and of the general mechanism of smoothing of singularities in 
string theory. We will relate the perturbative a' corrections to the integration over the first 
few excited string states (sect. 3.1); instead, the mechanism of singularity regularization is 
related to the exponential growth of the density of states; it is therefore basically due to very 
massive string modes. 

In sect. 4 we study the effect on the lowest order solution of the backreaction due to very 
massive string modes produced by the gravitational field. For technical reasons, that we will 
discuss in sect. 4, it turns out to be quite difficult to present an explicit non-singular solution. 
The result also suggests that treating massive modes production as a backreaction, in a full 
stringy regime, is not adequate. However, from the equations that include the creation of 
massive modes, we will be able to obtain at least a qualitative understanding of the role of 
massive modes production and of their possible regularizing effect. In sect. 5 we present the 
conclusions and we discuss possible phenomenological implications of a string phase with 
highly excited massive modes. 

2 Perturbative a' corrections 
2.1 The constant curvature solution 

Let us recall the results of ref. Including first order corrections in a', a possible form of 



the effective action, in the string frame, is [|I2|, [13 

S = / (f+^x^e 



(2.1) 
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where is the string length, ~ ^Jot' ■, k = 1, 1/2 for the bosonic and heterotic string, 
respectively, and we have neglected the antisymmetric tensor field. For type II strings 

T^ . Performing the field redefinition 



k 



0, and the first correction starts at order Rf,^p„ 



Qfiu giMu + ika')Sg^^, 0^0+ {ka')5(t) with 

= R^.u - 5^09.0 + ^7;..(V0)^ 5(t) = -\R+ {2d - 3)(V0) 



and truncating at first order in a' gives the action 
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(2.2) 



(2.3) 



where Rq^ 
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— 4i?^^ + R^ is the Gauss-Bonnet term. This form is particularly 
convenient because higher order derivatives in the a' correction cancel after integrations by 
part [|l^. We will discuss in some detail field redefinitions in sect. 2.2. 

We now restrict ourselves to a homogeneous and isotropic FEW background, ds'^ = 
N'^{t)dt^ — a?{t)d'x^ and homogeneous dilaton = 0(t). Varying the action with respect to 
a, 0, we get two dynamical equations of motion. Introducing H = a/a and specializing to 
d = ?) for notational simplicity, they read (in units ka' = 1): 



-6H{1 + H^) + 



1 + 



12H^ + met) 



- 12H{1 - H^) + 60(1 + H^) - ISH^ + 12H(p - 30^ + 
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(2.4) 



^ + 12H^<j)-6H^<j)'^ = 0. (2.5) 



The variation with respect to the lapse function gives, instead, a constraint on the initial 
values: 



6H'^ - 6H^ + 02 
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(2.6) 



The constraint is conserved by the dynamical equations of motion, and is therefore satisfied 
at any time if it is satisfied at the initial time. We can now try to solve eqs. ( |2.4| - |2.6| ) 
with the ansatz H = const = y, <p = const = x. The ansatz reduces the three differential 
equations to algebraic equations in x, y. At first sight we have three independent equations 
for two unknown variables. However, reparametrization invariance gives one relation between 
these equations such that if the constraint equation and, say, the equation obtained with 
a variation with respect to a (or, in the general anisotropic case, the d equations obtained 
with a variation with respect to ai, i = 1, ... ,d) are satisfied, then the equation obtained 
with a variation with respect to is automatically satisfied. 
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In the generic anisotropic case with scale factors ai [i = 1, . . . ,d), the ansatz therefore 
reduces the system of {d + 2) differential equations to {d + 1) algebraic equations in {d + 1) 
variables Hi = yi,<p = x, and this is true at all orders in a'. 

These algebraic equations are nothing but the requirement that there is a zero in the 
beta functional of the underlying sigma model, when the background is specialized to the 
form of our ansatz, and we will therefore write them as /3j(g) = 0, where i = 1, . . . , {d + 1) 
and g = {x,yi, . . .,yd). 

It is easy to generalize this result to the case of compact dimensions. Let us first consider 
the case in which the compact space is spatially curved, e.g. consider four- dimensional space 
times S"^, with metric 

3 

ds^ = N\t)dt^ - ^ a^{t)dx^, - b\t) {de^ + sin^ 9d(f)^) . (2.7) 
1=1 

(The extension to more general cases will be obvious.) Writing down the equations of 
motion we see that now b enters not only in the combination b/b, but also through terms 
~ l/b"^, which are due to the curvature of the sphere. Therefore the ansatz di/ai = constant, 
(p = constant can only be consistent if 6 = const, or b/b = 0. Again for the ansatz di/ai = 
yi,(j) = x,b = const the relation between equations of motion derived from reparametrization 
invariance eliminates one equation and we have a number of algebraic equations equal to the 
number of variables yi,x,b. 

The same happens if, instead, we compactify on a torus. In this case the explicit depen- 
dence on the scale factor, which forces it to be constant, comes from the winding modes, 
whose energy grows with the scale factor. 

In conclusion, the existence of a solution with Hi, (f) constant (and b constant, for compact 
dimensions) depends on whether the algebraic equations discussed above have real solutions. 
For the action (|2.3|) this is indeed the case, and the solution turns out to be an attractor 
of the lowest order pre-big-bang solution, which is therefore regularized 0. However, the 
inclusion of higher orders in a' or field redefinitions of the type used in eq. ( p^.2| ) produce 
different algebraic equations, which may or may not have real solutions. These issues, which 
were already noted in ref. 0, will be further discussed in the next section. 
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2.2 Scheme dependence of the results 

In principle, we would like to know the beta functions /?i(g) exactly, while what we have is a 
perturbative expansion in powers of a' . Somewhat optimistically, one might still hope to find 
a zero which, in units koi = 1, is located at gfj <^ 1, thus justifying a perturbative treatment. 
Unfortunately, an even more fundamental obstacle stands in the way. The problem is that, 
if we work at finite order in a', the coefficients of the algebraic equation, or of the expansion 
of the functions are subject to ambiguities. A straightforward way to understand this 
point is to observe that we can perform fields redefinitions that mix different orders in a'. 
The most general form of such redefinitions, at order a', is [0] g^^y g^y + {ka')6g^y, 
(j) ^ (j) + {ka')6(p, with 

Sg^,„ = aiR^y + a2di,(pdy(j) + a'ig^,y{V(t)f + a^g^^yR + a^g^^y^cj) , (2.8) 

5(l) = hiR + h2{S/(l)f + hU<p. (2.9) 
It is not necessary to include a term V^(9y0 in eq. ( p.8| ) because it can be reabsorbed by a 



general coordinate transformation [T^. After this redefinition, the new action, truncated at 
order a', is (in units ka' = 1) 

S = -J d^+^x^ e-<^ \r + (V0)2 - -^Rlyp^ + c^R^yR^" + csi?' + C3(V0)^+ 

+CiR'"'d^(t)dy(t) + C5i?(V(/))2 + c^RUcj) + c-jU(j){y(t)f + c^i^ct)?] , (2.10) 

with the coefficients q functions of aj,6j. (We have eliminated terms that can be reduced 
to the above terms after integrations by part or use of Bianchi identity). Despite the fact 
that we have eight coefficients q and eight parameters ai, . . . , 05, 6i, 62, ^s, we cannot fix the 
Ci at arbitrary values since, from the explicit expression of the q as functions of oi, . . . ,63, 
one finds that they satisfy a relation C2 + C3 + C7 + Cg = C5 + Cg. Within this constraint, 
however, the Cj can be chosen at will with the appropriate field redefinitions. The existence 
of a relation between the q means that there is a one-parameter family of field redefinitions 
which, at order a', leaves the action invariant. It is readily found to be 



= Cg,u {R - {Vct>f + 2Ucf) , (2.11) 

= C ( - ^(V0)^ + doA , (2.12) 



with ( a real parameter. We now ask whether the existence of a zero of the functions 
at a given order in a', is affected by the field redefinitions, and the answer is positive. For 
instance, we can fix ci = 1,C2 = — 1/4, C3 = 1/4, and vary C4 (which does not enter the 
relation between the Cj) setting all others c, to zero. For C4 = we have the action used 
in [0, and there is a zero of the beta functions. In = 3, for the isotropic case, it is located 
at X ~ 1.404,?/ ~ 0.616. Increasing C4 this zero disappears (escaping at infinity in the {x,y) 
plane) at a critical value C4 ~ 0.05. 

Therefore, while we expect that, at all orders in a', the existence of a zero in the beta 
functions is independent of field redefinitions, this is not true for the truncation at any finite 
order in a'. 

There are other useful ways to understand the existence of ambiguities in the perturbative 
expansion of the functions jSi, and these different points of view are believed to be equivalent. 
First, from the point of view of the underlying sigma model, they are due to the dependence 
of the perturbative coefficients of the beta functionals on the renormalization scheme. This 
dependence starts at two loops, i.e. from the terms ~ in the action, eq. ( 2.10 ). 

The effective low energy action, on the other hand, is constructed in such a way as to 
reproduce the string theory S'-matrix elements. From this point of view, the ambiguities 
in the coefficients come from the fact that some coefficients cancel in the computation of 
on-shell amplitudes [|T7i |T^ , p!8| . Suppose for instance that we want to fix the coefficients 
of the operators Rf^^p^,Rf^^, and R^ in the effective action. We would then expand these 
operators around the flat metric, g^^, = rj^y + h^^ and compute their contribution to the 3- 
point and 4-point amplitudes. However, the ~ part of ^—gR^^R^^ and of y/—gR^ vanish; 

therefore, the coefficients Ci, C2 in the action ( |2.10 ) cannot be determined from the knowledge 

2 



of the string amplitude with three on-shell gravitons. Both y/^R^^R^'^ and yJ—gR are 



non- vanishing when expanded at order /i^. However, the string amplitude includes one- 
particle reducible graphs with graviton exchange; to reproduce it we must therefore sum the 
contact and the exchange terms derived from the effective action. In the sum, the coefficients 
Ci,C2 cancel (see |jl8| for a detailed computation), and so they cannot be determined from 
the comparison with string amphtudes. The coefficient of -R^j^po- is instead fixed by this 
procedure, at the value —1/4. 
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2.3 Relation with scale factor duality 

Scheme dependence should not change physical results if we had the exact expression for 
the beta functions. It enters into play when we truncate at finite order in a'. When scheme 
dependence appears, some scheme will be 'better' than others, in the sense that the results 
obtained at finite order in this scheme will be closer to the exact result. Since scale factor 
duality plays an important role in motivating the cosmological model that we are discussing, 
one might hope that a scheme that respects scale factor duality at a given order in a' will 
be better, in the above sense. 

Scale factor duality has been generalized to O(a') in refs. [p!9| , pO| . The analysis of ref. |]19 
is very general, and refers to the full 0(d^ d) symmetry. The result of ref. [|19] is that there is 
one, and only one action invariant under scale factor duality at order a', and it is given by the 
action ( p.lQ|) with ci = 1,C2 = -1/4, C3 = 1/4, C4 = 1,C5 = -1/2, ce = 0, C7 = -1/2, cs = 0. 
At the same time, the duality transformation must acquire a' corrections. If we restrict to 
FRW metrics and vanishing antisymmetric tensor field, then ai 1/ai (or loga^ — loga^) 
must be generalized as 

log ai —(log a,,) — ka' ( — ) . (2.13) 

\aiJ 

The transformation of the dilaton is fixed from the condition that (p logc^i is invariant. 
In c? = 3, for isotropic metric, the equations f3i{x,y) = for this dual action have two 
real solutions]^ at (x ^ 1.526,?/ ~ 1.913) and at (x ~ 6.201,?/ ^ 1.931) (plus the solutions 
{—X, —y), which are always associated with the solution {x,y)). The existence of a pair of 
solutions (plus the sign-reversed pair) is related to the scale factor duality of the action. 
However, for the action that we are considering, scale factor duality is valid only at order 
a', and not exactly, since the transformation ( p. 131) , when applied to the 0{a') terms of the 
action, generates 0{a'^) terms that are truncated. Thus the two solutions are not exactly 
dual to each other. 

Integrating the equations of motion numerically, one finds that these solutions do not 
act as a late time attractor of the lowest order solution, which still diverges at some finite 
value of time. 

To further explore the relation with duality we have tried a different generalization of 



^Real solutions for this action exist for any d. This can be proved numericaUy up to large values of d, 
where the results matches with a large d expansion. 
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duality at 0{a'). In fact, the action found in |]T9l and the form of the duahty transforma- 
tion ( |2.13| ) are uniquely fixed only if we consider the general case of non-isotropic metrics. 
If however we restrict to isotropic metrics, a transformation 

loga ^ -(loga) - AfcaM H , (2-14) 

with arbitrary A, leaves the action ( p.lOj ) invariant at order a' if the coefficients Cj are Ci = 
1, C2 = -1/4, C3 = 1/4 + (A - 1)/18, C4 = 1 + (A - l)/9, c, = -cj2, cg = 0,cj = -2cs, Cg = 0. 
(The relation between the q is still satisfied.) For A = 1 we recover the action and the 
transformation of ref. [0. However, for every transformation ( p.l4| ) with given A we now 
have one action that is invariant .0 In a sense, this enlarged family of transformations is less 



interesting than the transformation found in [|T9l , because the latter is the only one that can 
be generalized to the anisotropic case and to the inclusion of the antisymmetric tensor field. 
However, this family of transformations also includes the case A = 0, i.e. the transformation 
a —>■ 1/a, without a' corrections, which is appealing for its simplicity. Furthermore, since 
this transformation does not generate higher orders in a', it is an exact invar iance of the 
action ( p.lO|) , with the appropriate choice of q, rather than an invariance at 0{a'). 

We have therefore studied the equations of motion for the action with A = duality. We 
have found that the functions Pi{x,y) have a pair of zeros (plus the sign reversed pair), and 
these zeros are related by exact duality invariance, as it should. But, again, these solutions 
do not act as late time attractors of the lowest order solution, which instead runs into a 
singularity. The type of singularity is however different from the one encountered for the 
case A = 1. It is interesting to understand in some detail what happens. The dynamical 
equations of motion can be written in matrix form 

-'(:)■(;:) 

where A is a 2 x 2 matrix that depends on H, 0, and /i, /2 are functions of H, 0. At a finite 
value of time we find that det A = 0, and therefore the equations become singular. As a 
simpler example of a similar situation, consider the differential equation {1 — h)h = h for 
•^Basically, we have a much less rigid structure because many expressions which are independent in the 
anisotropic case, such as J^i ^"^^ (X^i , cohapse to the same expression in the isotropic case, and we 
can arrange cancellations between them. 
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some function h{t). For small h the solution grows like e*, until the 'determinant' (1 — h) 
vanishes, so the situation is similar to our case. Here however we can integrate the equation 
exactly, and find t{h) — log/i — h. Above a critical value t — —1, this relation cannot 
be inverted while, below this value, there are two branches h+{t) and h^{t) that merge at 

This suggests that what happens in our case is that the doubling of solutions due to 
duality transformations with A = produces pairs of solutions of the dynamical equations 
of motion, which merge at some critical value of time, and thereafter move into the complex 
plane. (These solutions are not dual to each other. One corresponds to a small initial value 
of H and H > 0, and the other to H large, and H < 0, while duahty changes H —H). 
If this conjecture is correct, then it appears that duality works against the regularization 
mechanism. In any case, neither a scheme that respects duality with A = 1 nor a scheme 
that respects duality with A = provides a realization of the regularization mechanism, at 
0{a'). 

3 The physics behind the regularization of the singu- 
larity 

The conclusion of the previous sections is that perturbative a' corrections can in principle 
regularize the singularity; however, to determine whether this actually happens we should 
know the functions Pi at all orders in a', because results at finite order are scheme dependent. 
The attempt to fix the scheme using scale factor duality as a guiding principle does not 
give encouraging results. We therefore ask whether there are general physical principles 
that suggest that the singularity is indeed regularized. Generally speaking, we expect that 
string theory eliminates all unwanted singularities. However, is the physical mechanism that 
ehminates the singularity related to perturbative a' corrections? Or should we look for a 
different type of corrections? In this section we discuss this question. 

3.1 Effective action and massive string modes 

As it is well known, the low energy action of string theory can be obtained either by com- 
puting the beta functions of the sigma model or from the requirement of reproducing the 
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S-matrix elements. The two methods are equivalent, as was proved in Still, the phys- 
ical pictures underlying the two computations are quite different. In particular, in the first 
method the massive string modes seem to play no role at all. 

In fact, what one does in order to include massive modes at the sigma model level is the 
following. The sigma model action is written including the interaction with the background 



fields representing the massive modes P^ , 
S =^J rfVv^ 



+-K,,^{x)h'^h^'d,x^'djX''dkXPdiX^ + . . 



a 



(3.1) 



where together with the metric, dilaton and antisymmetric field, we have displayed an ex- 
ample of a field at the next mass level. Once we include massive fields at the first excited 
level, the whole tower of massive fields must be included, since they are associated with 
non-renormalizable interactions from the point of view of the 2-dimensional theory, and all 
possible terms are in principle generated by the sigma model loop (i.e. a') expansion. 

However, when we compute higher loops in the beta functionals of the massless modes, 
the massive modes play no role: the operators associated with massive modes are non- 
renormalizable and do not mix with the massless sector. More generally, to renormalize 
background fields at a given mass level, we need only fields at the same or lower mass 



level [23 



However, since in string theory we have an infinite tower of massive modes, while the 
effective low energy action is written in terms of massless fields only, it is clear that we have 
somehow integrated over the massive modes. This physical interpretation is confirmed by 
the computation of the effective action through the comparison with S-matrix elements. In 
this case the effective action, at all orders in a' but without e*^ corrections, is obtained by 
requiring that it reproduces the string amplitudes at genus zero, with an arbitrary number of 
insertions of vertex operators for the massless fields. In the field theory language, we expect 
that the a' corrections to the amplitudes can be represented as a sum of tree level exchange 
graphs with massive string modes in the intermediate state. This is indeed the case, as can 
be read from the computation of Gross and Witten [|T^ of graviton scattering in type II 



superstring theory. In this case the tree amplitude is ||2^, |T4 

r(-s/8)r(-t/8)r(-M/8) 



A 



N 



128 



r(i + |s)r(i + |t)r(i + |«) 
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(3.2) 



where is the gravitational constant, s, t, u the Mandelstam variables and K are kine- 
matical factors. 

Expanding the term in brackets for small momenta, one gets a leading term ~2^/{stu), 
which, when inserted back in eq. (p^), gives the tree level scattering amplitude derived 
from the lowest order effective action; and a correction term —2({3), which gives the lead- 
ing correction to the effective action; from the kinematical factors K, K one finds that it 
corresponds to a term in the effective action quartic in the Riemann tensor (remember 



that for type II superstrings the coefficient k in eq. (|2.1|) is zero, i.e. there is no ~ R^^upa 
correction). 

Writing the Riemann zeta function as C(3) = Y^^=if^~^i we see that the correction can 
be interpreted as a sum over an infinite tower of intermediate states, and that states at level 
n give a contribution ~ 

This means that the a' corrections are dominated by the first few massive string states. 
For instance, summing over the first 10 excited states, we get E^ii 1/^^ - 1-1975 . . ., to be 
compared with the value of the Riemann zeta function (^(3) ~ 1.2020 . . .. 

Having obtained a physical picture of the mechanism that is responsible for the pertur- 
bative a' corrections to the effective action (integration over the first few massive modes), 
we proceed in the next section to discuss the physical mechanism that is responsible for the 
regularization of the singularity. 

3.2 Small distances vs. large energy singularities 

It can be useful to distinguish between small distance singularities and large energy singu- 
larities. In field theory small distances means large energies, and the distinction is therefore 
meaningless. In string theory, however, the behavior of amplitudes in high energy scattering 



suggests the existence of a generalized uncertainty principle of the form p5 



Ax > + const. GivAp. (3.3) 

(Such an uncertainty principle is also suggested by general properties of quantum black 
holes |2^.) Therefore in string theory the connection between high energies and small 
distances is not completely trivial and it is useful to consider the two cases separately. 

The most basic reason why we expect strings to regulate singularities at small distances 
is that the very notion of invariant point-like event is not meaningful in string theory. While 
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in quantum field theory an invariant event can be defined through the sphtting of a particle 
in two, no such concept exist for strings, since the point in space-time where a string splits 
in two depends on the Lorentz frame used p7 |. 

At large energies, the regulating mechanism is instead the exponential growth of the 
density of states. The asymptotic density of states has the form 

rf(M)~ f— V'e*^/^% (3.4) 



VMo/ 

where Mq = cj y/a' and 6, c are numerical constants, which depend on the string theory under 
consideration. An example of how the density of states prevents quantities with dimension of 
energy from growing arbitrarily large is given by the fact that we cannot raise the temperature 
of a system beyond the Hagedorn temperature, Thag = because otherwise the canonical 
partition function diverges. 

Another interesting example is given by the existence of a limiting value of the electric 
field for an open bosonic string. In the presence of an external electric field the rate for 



charged-string pair production is in fact p8 



E9^E(-i)Mx ^-p -^(^1 + ^') • (3-5) 



w 

S k=l 



The sum goes over all physical string states S, with mass M5, and qs is a factor that depends 
on the charge of the state S\ in the weak field limit, e ~ eE + o(£'^), where e is the total 
electric charge of the string, and one recovers the Schwinger result. As long as e is finite, the 
factor exp(— 7rA;M|/|e|) ensures the convergence of the rate, even if the sum over the states 
5, for large masses, is an integral with the exponentially growing density of states (|3.4|) . 
However, at a critical value of the electric field, e goes to infinity and the rate therefore 
diverges. 

The cosmological singularity is a large energy singularity: the Hubble constant, and 
therefore the curvature or the energy density, diverges on the lowest order solutions of the 
equations of motion. We therefore expect that the regulating mechanism has to do with the 
existence of an infinite set of massive states; in fact, we expect something similar to what 
happens in eq. ( p.5|) , since the existence of a maximum electric field and of a maximum 
gravitational field do not seem to be two fundamentally different problems. 

So, the first conclusion is that we should not necessarily expect that perturbative a' 
corrections regulate the singularity. As we discussed in sect. 3.1, these corrections can be 
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accounted for with an accuracy better than one per cent summing over the first ten mass 
levels. The regularization of the singularity, instead, must be crucially related to the existence 
of an infinite tower of massive states. Of course, this does not mean that the mechanism 
proposed in cannot regulate the singularity. It simply means that we do not have any a 
priori argument that allows us to confidently state that the mechanism will operate. 

The second point suggested by the above considerations is to look for massive modes 
production in a gravitational field. This will be discussed in the next section. 



4 Production of highly excited modes 

A time- varying gravitational field produces particles. The general mechanism is the existence 
of a non-trivial Bogoliubov transformation between the in and the out vacuum. For instance, 
in De Sitter space with Hubble constant H, a detector moving on a geodetic sees a thermal 
bath of particles at a temperature 

T{H) = ^. (4.1) 

This suggests that in De Sitter space a massive string mode with mass M is created with 
a probability proportional to exp{—M/T{H)). One might wonder whether a truly stringy 
computation of massive modes creation could give a different answer from the field theory 
result ([4.1|). However, Lawrence and Martinec p9[ (see also [^) used string field theory to 



compute massive modes production in a FRW space with 

ds^ = C{r]) (drf - dx^) , C(r/) = A + B twh{pr]) , (4.2) 
and they obtained the same result that holds in this model in the field theory limit, discussed 



e.g. in ref. |^| 



The above result suggests that the rate of energy density production due to massive 
modes creation during the expansion is approximately 



P 



1 

j^J^^ dM Md{M)e-^^^^/^ 

dxx^-''e-^% (4.3) 



s 



where 

,.M^_1.JI£,_1 (4.4, 
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and d{M) is the density of states, given asymptotically by eq. (|3.4| ). This formula, which has 
been suggested in [^] as valid for all the mass levels, is also qualitatively confirmed by the 
explicit computation of the Bogoliubov coefficients for the graviton production, i.e. at mass 
level zero [0, 32]. In this case, in fact, one finds that is a number of order 1 for energies 
of the order of the string mass (red-shifthing the frequency at the present epoch, this means 
that = 0(1) if the physical frequency, as seen today, is of order 10 GHz), and has an 
exponential cutoff for larger frequencies. Actually, the spectrum also shows some deviations 
from an exact black body spectrum (it goes like for low frequencies but as u^^'^'^, where 
/i depends on details of the string phase, for intermediate frequencies.) These deviations, 
however, are not very relevant to the present discussion.0 

Equation ( ^.3| ) clearly has a regularizing effect on the growth of H. When H exceeds a 
critical value He, given by He/ (27r) = Mq, the production rate diverges. This result can also 
be interpreted by noting that H/ (27r) is the Hawking temperature in De Sitter space, and 
therefore the condition on H^ is the requirement that the Hawking temperature does not 
exceed the Hagedorn temperature, Thag = Mq. 

The effect that we are discussing is non-perturbative in a'. In fact, for 7 > the integral 
in eq. ( [4.3|) is an incomplete gamma function, and from its known asymptotic expansion we 
find that, in the limit a' ^ 0, 

,.(Lf\.J-4^], (4.5) 



.a' J [ HVa' ^ 

and it is therefore non-analytic in a' at a' = 0. 

To understand in more detail the evolution of the cosmological model, we tentatively 
use the expression ( ^.3| ) for the energy density produced, as a backreaction term in the 
cosmological equations. The equations of motion, without perturbative a' corrections, and 
including matter sources with = {p(t), —6jp{t)), are 0, |] (in the isotropic case, in d 
spatial dimensions) 

0^ - dH^ = 2Af-^e^p (4.6) 
20 - 0^ - dH'^ = (4.7) 



^In an expanding Universe one should also add to p a term —aHp representing the dilution of string 
density due to the expansion [^3|, and a = d — 1 for long strings. However, this term is not very important 
for the qualitative discussion presented below, as we checked also numerically. 
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H-H(l) = \i-^e^p (4.8) 
where (p = (p — d log a, e'^p = e'^p, e'^p = e'^p. We use eq. ([4 .31) for the time derivative of p; the 



pressure p is determined by the conservation of the energy momentum tensor, which gives 

p + dH{p + p) = 0. (4.9) 

This is equivalent to requiring that the constraint equation ( [4.6| ) is conserved by the dy- 
namical equations of motion ( [4.7|) and ( [4.8|) . Since we know p rather than p,p, these are in 



principle integro-differential equations. However, we can combine eq. ( [4.8|) and eq. ( [4.6|) so 
that only the combination p + p appears, which is expressed through p using eq. ([4.9|) , and 
we are left with ordinary differential equations. Writing everything in terms of (p rather than 
0, the dynamical equations of motion can then be written as 

2if + 02 - 2{d + 1)H<P + d{d + l)H^ = --X'^e^'^ , (4.10) 

d H 

20 - 2dH - 0^ + 2dH^ - d{d + l)H'^ = . (4.11) 

The constraint equation 

+ d{d-l)H^ -2dH^ = 2Xi'^e'^p (4.12) 

is still an integro-differential equation. However, we impose the constraint at the initial 
time to ^ ~1 (in units = 1), when the right-hand side is negligible. The constraint is 
then automatically preserved by the evolution (we use its conservation as a check of the 
numerical integration). We have then integrated these equations numerically. The results 



are qualitatively the same for any d and any value of 6, c in eqs. ( [4.3|) and (4^) that we have 
tried. The plots that we present refer to d = 3 and to the values 6 = 10, c = 10 (this value 
of c is responsible for the values ~ 0(10) of the H-scale in Fig. 1). Similar results holds for 
c = [(2 + \^)n]^^ ~ 0.09, which is the value for the heterotic string. 

Figure 1 shows the behavior of the Hubble parameter versus cosmic time t. At large 
negative values of t we recover of course the lowest order solution, since the production of 
massive string states is exponentially small; H{t) then reaches a maximum value and then 
very abruptly bounces off and goes to zero at a finite value of t. At this point the numerical 
integration stops.0 The behavior of is shown in fig. 2, while fig. 3 shows p{t). 

^The integration stops because the integral in eq. (^^) is only defined for 7 > 0. For negative values of 

15 



These plots show that some important aspects of the regularization mechanism are still 
missing from our equations. In fact, even if H does not diverge, H becomes very large when 
H goes to zero, and the Ricci scalar TZ = —6{H + 2H^) in this regime is dominated by H and 
is large. The missing ingredient can be traced back to eq. ( |4.1| ). This form of the temperature 
holds for De Sitter space, where H is constant and H = 0. It is inappropriate at a stage of 
the cosmological evolution when H is small compared to the critical value but H is large. 
Indeed, at this stage our equations predict no massive mode creation, while a large H, as any 
strong time-varying gravitational field, will produce particles, much the same as a large H. 
This is quite clear physically, but it is anyway interesting to check this point in a toy model. 
One can consider for instance a model with a scale factor such that H = A/ cosh (cut), with 
A,uj > 0. If y4 <^ 1 but Auj ^ 1 the Hubble constant is always small during the evolution 
while H reaches a maximum value -ffmax = Aij/^/2 ^1. A straightforward computation 
gives a production rate of massive particles ~ exp(— M/T) with T ~ i/max- 

So, at least at the level of phenomenological description of the string phase, one should 
generalize eq. (|4.1|) to 

T = T{H,H,H,...,^,4>,...), (4.13) 

where all higher order derivatives of H appear. Note that we also included a dependence on 
<p because, changing frame, say from the string to the Einstein frame, the derivative of the 
dilaton field contributes to the gravitational field, and therefore in a generic frame even <p 
must contribute to massive modes production. For physical reasons, we expect that T will 
be a sum of positive definite contributions. For instance, a possible plausible form could be 



T = — iH' + ciH' + C2r + ...) , (4.14) 

with ci,C2 > 0. Now it becomes technically impossible to solve the differential equations, 
since already including a term H into T, the highest derivative of H in the differential equa- 
tions appears as argument of the incomplete gamma functions, which makes the equations 
intractable, even numerically. However, the qualitative effects of eq. ( |4.1(J| ), combined with 



an expression of the form of eq. ( |4.13|J4.1^ ) for the temperature, can be understood as follows. 



H we could replace H is eq. (4.4) with \H\. This is motivated by the fact that, for instance, in the model 



defined in eq. (4.2) the production rate is the same independently of whether we consider a model with a 
scale factor increasing in time or decreasing in time with the same rate. However, as we will now discuss, a 
more crucial modification of our equations is needed. 
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When any of the derivative oi H oi (j) exceeds a critical value of order one in units = 1, 
the right-hand side of eq. ( |4.1U| ) diverges, while the left-hand side is finite. So there can be 



no solution where H, H,(j), . . . are larger than their critical value. Similarly to what happens 
to H in fig. 1, when H, 0, etc. approach their respective critical values, the effect of the term 
p in eq. ([4.10|) is to repel them away. Furthermore, H cannot go to zero. Indeed, unless all 
higher derivatives H,(j), . . . vanish (as it happens at t ^ — oo if the evolution starts from 
the string perturbative vacuum), then p{H, H,(f), . . .) does not vanish for H —>■ 0, contrarily 



to what happens when we use the expression T{H) = H/27r, see eq. ( |4.5| ). Therefore the 
right-hand side of eq. ( [4.10| ) diverges if — > while the left-hand side is finite, because 
H, <p are bounded. Therefore, there can be no solution with H approaching zero. So, we 
expect a solution in which the scale factor in the string frame is always expanding, H > 0, 
possibly in a very complicated way, but still without divergences in the curvature or in any 
other quantity.0 

The results obtained trying to incorporate massive mode production as a backreaction 
on the metric and dilaton field show that massive modes may have a very important effect 
on the evolution of the cosmological model. However, the very complexity of the equations 
that describe the backreaction of massive modes on the massless modes 0, g^i, indicates 
that a treatement as backreaction is not appropriate. If we enter a full stringy regime, 
with an infinite tower of excited massive fields, it is not surprising that a description of the 
cosmological model in terms of classical evolution of massless fields runs into difficulties. 

Finally we stress that, even if our results allow to obtain some understanding of the 
physical mechanism that regularizes the singularity of the lowest order solution, still this 
mechanism cannot be the full story, since it does not provide a solution that, besides be- 
ing possibly non-singular, is also smoothly connected with the present radiation dominated 



phase. This is the so-called graceful exit problem [H, pi ITO|, ITl|. It can be posed as follows. 



From eq. ( [4.6| ) one sees that there are two solutions for 0, given by 

= ± [dH^ + 2\i~^e^p) ^'^ . (4.15) 

These two solutions are usually referred to as the (-I-) and (— ) branches, depending on the 
®We stress that this is not a proof that there will be a non-singular solution. In principle it may happen 
that at some value of t two different solutions 'collide' and move into the complex plane, as we have seen in 
sect. 2.3 for the equation with perturbative a' corrections. 
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choice of sign in front of the square root. The pre-big-bang solution starts on the (+) branch; 
this ensures that the solution start from a regime of low curvature and couplings, and goes 
through a phase of inflation driven by the kinetic energy of the dilaton field. The final state, 
however, should correspond to a radiation dominated FRW model with constant dilaton, 
and therefore = — dH < 0, and we are on the (— ) branch. Therefore at least one branch 
change must take place. Eq. ([4.15|) shows that a necessary condition for this to happen is 
that, at some point in the evolution, the matter density p, in the string frame, be negative. 
Further global conditions are necessary to ensure that no further branch change back to the 
(+) branch takes place, and have been discussed in refs. |1T0| , p!T[] . However, we see that, if we 
include massive modes as the only source of energy density, p is necessarily positive and no 
branch change can occur. Other source of corrections, like a' correction, provide an effective 
contribution to p which is negative and therefore allows branch changes. Still, to satisfy 
the global conditions for a succesful exit completion, e*^ corrections must also be included. 
An example in which the combined effect of a' corrections and corrections provides a 



successful exit has been discussed in [|rT|. So, it is clear that the effect that we have studied 
in this paper is only an ingredient of the mechanism that finally can produce a well-motivated 
and non-singular cosmological model, that eventually matches standard cosmology. 



5 Conclusions 

The lowest-order cosmological solutions obtained from the string effective action are singular, 
and therefore corrections to the lowest order action must play a crucial role. In string 
theory there are two kind of corrections, 'stringy corrections', parametrized by a', and loop 
corrections, parametrized by e*^. In this paper we have drawn attention to the fact that a' 
corrections come into two types, perturbative and non-perturbative, and that they have a 
distinct physical origin. Perturbative corrections are related to the integration over the first 
few massive modes while non-perturbative corrections comes from the exponentially growing 
density of massive states. As we discussed in sect. 3, this physical interpretation makes clear 
that non-perturbative a' corrections are expected to have a regularizing effect. 

The explicit verification of the regularization mechanism meets difficulties both for per- 
tubative and non-perturbative a' corrections. In the former case, a truncation at any finite 
order in perturbation theory gives scheme-dependent results, so that an all-order computa- 
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tion is in principle required; since an all-order perturbative computation is anyway faced with 
the problem of the convergence and resummation of the asymptotic series, this is already a 
hint of the fact that the full answer requires the inclusion of non-perturbative effects. 

In the case of non-perturbative corrections related to massive modes production, the 
technical difficulty that we have to face is the inclusion of the effect of all higher order 
derivative of H and in the computation of the production rate of massive modes. However, 
independently of the technical difficulties, the mechanism has a clear physical interpretation, 
which suggests that it regularizes the lowest order solution. 

We conclude by pointing out briefly a possible phenomenological consequence of the 
scenario in which in the string phase we have highly excited string modes. The massive 
string states produced will eventually decay. The decay of the excited string levels has been 
studied in ref. |3|]. The result is that, for large level number N , the width is dominated by 
decays in which one of the products is massless, i.e. by transitions of the form N N — 1 
with emission of a graviton or another massless particle. From the mass formula Mjj ~ A^, 
the level spacing is AMn ~ I/Mn- If a string at an excited level decays at level — 1 
with emission of a graviton of frequency u = AMn, we get transitions at the frequencies 
In = fo/VN. The most naive estimate of the frequency /o, as seen today, is obtained by red- 
shifting the string mass from a string time to the present epoch, which gives /o — '-'^today/ (^Tt) 
of the order of 10 — 10^ GHz, which is also the typical cutoff value of the gravitons produced 
by quantum vacuum fluctuations 0, However, this radiation is produced within the 

string phase, rather than at the transition between the string phase and the FRW phase, 
and a further red-shift must be included, to take into account the expansion from the time 
of production to the onset of a standard FRW cosmology. We cannot compute it, in the 
absence of a knowledge of the mechanism that terminates the string phase. There is, however, 
the interesting possibility that, because of this further redshift and/or the factor in 
/at, this gravitational radiation might be shifted in the frequency band accessible to the 
LIGO/ Virgo gravitational wave detectors, which are expected to operate in the 6 Hz-1 kHz 
region. 

Concerning the intensity of this gravitational radiation, it can be estimated as follows. 
If there is not a subsequenty inflationary phase below the string phase, also the photons 
that we observe today in the CMBR have been produced during the large curvature phase 
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of the model, by massive modes decay. Since at this scale the rates for photon and graviton 
production are comparable, this suggests that the overall intensity of the relic gravitational 
wave background produced by the decay of massive string modes has an energy density 
comparable to the energy density of the 2.7 K radiation. 
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Figure 1: The Hubble parameter H as a function of time. The initial conditions for the 
numerical integration are to = — 20, 0(to) = —10, H{to) = — l/(to\/3); (j){to) is determined 
by the constraint equation. 
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Figure 2: as a function of time. Same initial conditions as in fig. 1. The integration stops 
when H reaches zero. 
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Figure 3: p as a function of time. Same initial conditions as in fig. 1. 
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